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ABSTRACT 

This paper presents a simple unknown input estimation technique using generalized matrix inverse. The full 
order observer constructed by g-inverse is extended and implemented for this purpose. The necessary and 
sufficient condition has been introduced for existence of such an estimator. The proposed method is illustrated 
by numerical example (Two loop autopilot in pitch plane) and the MATLAB simulation results. 
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I. INTRODUCTION 

Unknown input estimation has a special importance in many industrial applications like fault detection 
and identification, cryptography or parameter identifiability, track seeking, railway signaling etc. Several 
researchers have designed different observers that have been extended to develop UIO and to estimate unknown 
input. D. G. Luenberger first designed an observer in 1971 which approximately reconstructed missing state 
variable information [3]. Prof. Das and Prof. T. K. Ghoshal jointly proposed the method of construction of 
reduced order observer by using generalized matrix inverse [4]. Now to estimate the unknown input, it is 
necessary to design an unknown input observer (UIO). In [5], a constructive solution to the problem of 
designing a reduced order Luenberger observer for linear systems subject to arbitrary unknown inputs has been 
presented. A direct full-order observer design procedure for a linear system with unknown inputs is presented in 
[6]. Avijit Banerjee, Partha Pratim Mondal and Prof. Gourhari Das constructed a full order observer using the 
concept of generalized matrix inverse [10]. Alexander Stotsky and Ilya Kolmanovsky presented a paper on 
unknown input estimation techniques for automotive applications [11]. Unknown input estimation technique for 
linear discrete-time systems has been given in [7]. The problem of simultaneous estimation of states and 
unknown inputs for a class of non linear chaotic systems has been addressed in [12]. Avijit Banerjee and Prof 
Gourhari das presented a method of estimating the unknown input using reduced order Das and Ghoshal 
observer [4] in [13]. However, in this paper the full order observer constructed in [10] has been extended for 
estimation of unknown input. 

Notations: In this paper, R n and R nxm denote the n dimensional Euclidean space and the set of all nxm real 
matrices respectively. / and are the identity and null matrices respectively with appropriate dimensions. The 
superscripts "T" represents the transpose of a matrix and the superscript "g" indicates the Moore-Penrose 
generalized matrix inverse of a matrix [8], [1]. 

II. MATHEMATICAL PRELIMINARIES 



If A £ R mxn is a matrix then there exists a unique matrix A 9 £ R nxm , which satisfies the following 
conditions: 

(AA 8 ) T = AA 8 (1) 
(A g A) T = A E A (2) 
A A g A = A (3) 
A g A A E = A s (4) 

Consider a system described by linear equation, 

Ax=y " (5) 

where matrix A £ R mxn , known vector y £ R m and unknown vector x £ R n . Eq. (5) is consistent if and only if, 
AASy =y (6) 

Now, if eq. (5) is consistent then the general solution of eq. (5) is given by 

x = A3y + (I - A 9 A)r (7) 
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([1] Graybill 1969 p. 104). Where r £ W is an arbitrary vector. 

III. PROBLEM FORMULATION 

Consider an LTI system described by 

x = Ax + Bu + Ev (8) 
y = Cx (9) 

where the state vector x £ R n , known input vector u £ R ml and the unknown input vector v £ M. m2 .y £ R m is 

the output vector. A, B, C, and E are known constant matrices with appropriate dimensions. It is assumed that 

the system is state observable. 

The system dynamics in the absence of unknown input can be represented as 

x u = Ax u + Bu (10) 
y u = Cx u (11) 

Now, subtracting eq. (10) from eq. (8) and eq. (11) from eq. (9) respectively the state space description of the 

auxiliary system can be obtained as, 

x v = Ax v + Ev (12) 
y v = Cx v (13) 

where x v — x — x u is the system response due to unknown input only and y v — y — y u is the output response of 

the system due to unknown input only. 

The state of the auxiliary equation can be obtained from eq. (12) and (13). Now the general solution of eq. (13) 
is 

x v = C 9 y v + Q -C 9 C)h (14) 
where h £ R n whose elements are arbitrary function of time. 
Form eq. (12) and (13) we have 

y v = CAx v + CEv (15) 
From eq. (14) and (15) we have 

y v = CA(I - C 9 C)h + CAC 9 y v + CEv (16) 
Putting x v from eq. (14) into eq. (12) we get, 

(/ - CSQh =AQ- C 9 C)h + AC 9 y v + Ev - C 9 y v (17) 
The general solution of eq. (17) for h is given by, 

h = (I - C 9 C)A(l - C 9 C)h + Q - C 9 C)AC 9 y v + Q- C 9 C)Ev + C 9 Cz (18) 
where z £ R n is a vector whose elements are arbitrary function of time. 
From eq. (18) and (16), the observer dynamic equation can be written as, 

ti = [(I - C 9 C)AQ - C 9 C) - KCAQ - C 9 C)]h + [Q - C 9 C)AC 9 - KCAC 9 ]y v + [Q - C 9 C)E - KCE]v + 

C 9 Cz + Ky v (19) 
where K is the observer gain and h -* h as x -> x. 

In order to eliminate the first derivative of y v explicitly present in (15) the following substitution has been made. 

h = q + Ky v (20) 
and the observer dynamics can be written as 

§ = {(/- C 9 C)AQ -C 9 C)- KCAQ - C 9 C)}q + {(/ - C 9 C)E - KCE}v + 

1(1 - C 9 C)AC 9 - KCAC 9 + Q - C 8 C)AQ - C 9 C)K - KCAQ - C 9 C)K}y v + C 9 C (21) 
Without loss of generality, the arbitrary vector z may be chosen as a null vector for simplicity and less 
computation, the observer dynamics can be represented as, 

§={(/- C 9 C)AQ - C 9 C) - KCAQ - C 9 C)}q + {(/ - C 9 C)E - KCE}v + 

{(/ - C 9 C)AC 9 - KCAC 9 + Q - C 9 C)AQ - C 9 C)K - KCAQ - C 9 C)K}y v (22) 

Condition for existence of UIO: 

To nullify the effect of the unknown input from observer dynamics the observer gain parameter K can be 
designed such that 

(/ - C 9 C)E - KCE = (23) 
The general solution of K is given by 

K= Q -C 9 C)E(CE) 9 +K V {I -(CE)(CE) 9 } (24) 
where K v is an arbitrary matrix and the arbitrariness of K is also depends on K v . Putting the value of K from eq. 
(24) into eq. (23) we get, 

(/ - C 9 C)E(CE) 9 (CE) = Q - C 9 C)E (25) 
Eq. (25) is the consistency condition of eq. (23). Putting eq. (24) into eq. (22), it can be shown that the observer 
matrix would be (A,, - K V C V ) where, 

A v = (7 - C 9 C)A(I - C 9 C) - (I - C 9 C)E(CE) 9 CA(I - C 9 C) (26) 
and C v = {I - (CE)(CE) 9 }CA(I - C 9 C) (27) 
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The arbitrary matrix K v should be chosen such that the real part of the eigen values of observer matrix 
{A v — K V C V ) become negative. 

Then the UIO for auxiliary system, given by eq. (12) and (13) can be obtained as 

3 = (A v - K v C v )q + {(/ - C 9 C)AC 9 - KCAC 9 + (/ - C 9 C)AQ - C 9 C)K - KCAQ - C 9 C)K}y v (28) 
and the estimated state of the auxiliary system can be expressed as 

% v = Q-C9Qq + {C9 + (J-C9C)K}y v (29) 
The general solution for estimated unknown input can be expressed from eq. (12), 

v = E 9 (Jc v - Ax v ) + (/ - E 9 E)h (30) 
where h is any arbitrary vector. From the theory of generalized matrix inverse it can be concluded that if the 
matrix E is of full rank the unknown input can be estimated as 

v = E 9 (% V -Ax v ), since E 9 E = /. (31) 



IV. NUMERICAL EXAMPLE 

Consider the system of a class of flight path rate demand missile autopilot in state space model as 
described in [9]. In this example the states of the state vector x are: 
x x = Flight path rate(y) 
x 2 = Body rate(g) 
x 3 = Elevator deflection(r|) 
x 4 = Rate of change of elevator deflection(?7). 

The state space model of a classical two loop autopilot has been given below: 
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The following numerical data for a class of guided missile have been taken for MATLAB simulation: 
7>2.85 sec; a 2 =0.00142 sec 2 ; o h =5.6 rad/sec; ( a =0.6; K b = -0.1437 per sec; v = 3000 m/sec; K p = 28.99; 
K q = -1.40; <o a =180 rad/sec; 



We have taken the unknown input matrix, E = 



and the arbitrary matrix K v is chosen as, k v — 



The above system is observable and the existence condition of unknown input has been satisfied. 
For simulation the known input is taken unit step and the unknown input is taken as v — 500e _50t stnl000t. 
Simulation responses for estimation of states and the responses of the system with and without unknown input 
are given in Fig. 4. 1-4.4 below, where black doted lines denote the estimated signals, red firm lines indicate the 
response of the state with unknown input and green firm lines indicate the response of the state without 
unknown input. 

In Fig. 4.1 to 4.4, x, denote the system states with unknown input, xu, denote system states without 
unknown input and x t denote corresponding estimated states with unknown input for i =1, 2, 3, 4. In Fig. 4.5, v 
denotes unknown input and v denotes corresponding estimated signal. From the simulation responses it can be 
seen that the estimated signal tracks well the high frequency unknown signal ' v', i.e. 'D' follows V. 
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Fig. 4.1: Flight path rate 
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Fig.4.2: Body rate 
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Fig. 4.3: Elevator deflection 
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Fig. 4.4: Rate of Elevator deflection 
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Fig. 4.5: Unknown input 



V. CONCLUSION 

The unknown input estimator has been presented for LTI system extending the full order observer as 
reported in [10]. In this process the system is firstly transformed to an auxiliary system. Then using standard 
methodology an observer has been developed for the auxiliary system. Now by applying simple algebraic 
equation, estimated unknown input is obtained. The necessary conditions are proposed and solved using 
generalized matrix inverse. Illustrated numerical example with simulated results show the estimated unknown 
input which tracks the actual unknown input well enough, even if the signal frequency is very high. 
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